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We study graphene which has both spin-orbit coupling (SOC), taken to be of the Kane-Mele form, 
and a Zeeman field induced due to proximity to a ferromagnetic material. We show that a zigzag 
interface of graphene having SOC with its pristine counterpart hosts robust chiral edge modes in 
spite of the gapless nature of the pristine graphene; such modes do not occur for armchair interfaces. 

Next we study the change in the local density of states (LDOS) due to the presence of an impurity 
in graphene with SOC and Zeeman field, and demonstrate that the Fourier transform of the LDOS 
close to the Dirac points can act as a measure of the strength of the spin-orbit coupling; in addition, 
for a specific distribution of impurity atoms, the LDOS is controlled by a destructive interference 
effect of graphene electrons which is a direct consequence of their Dirac nature. Finally, we study 
transport across junctions which separates spin-orbit coupled graphene with Kane-Mele and Rashba 
terms from pristine graphene both in the presence and absence of a Zeeman field. We demonstrate 
that such junctions are generally spin active, namely, they can rotate the spin so that an incident 
electron which is spin polarized along some direction has a finite probability of being transmitted 
with the opposite spin. This leads to a finite, electrically controllable, spin current in such graphene 
junctions. We discuss possible experiments which can probe our theoretical predictions. 

PACS numbers: 73.20.-r, 73.40.-c, 73.63.-b 


I. INTRODUCTION 

The last several years have witnessed a tremendous 
amount of research on graphene, both theoretical and 
experimentalip— . Graphene is a two-dimensional hexago¬ 
nal lattice of carbon atoms in which the tt electrons hop 
between nearest neighbors. At half-filling, the spectrum 
is gapless at two points (called K and K') in the Bril- 
louin zone, and the energy-momentum dispersion around 
both those points has the Dirac form Ej: = hv\k\, where 
V ~ 10®to/s is the Fermi velocity. The Dirac nature 
of the electrons gives rise to many interesting properties 
of this material, such as Klein tunneling through a bar¬ 
rier—, novel effects of crossed electric and magnetic fields^, 
qualitatively different transport characteristics of super¬ 
conducting graphene junctions^^— , possibility of multi¬ 
channel Kondo physicsi^rJ^, interesting power laws in 
the local density of states (LDOS) induced by an impu- 
rityil^^, and atomic collapse in the presence of charged 
impurities2i. 

Recent years have also seen extensive research on topo¬ 
logical system o^^i^^ . These systems have a bulk spec¬ 
trum which is gapped; however, the topological prop¬ 
erties of the bulk states ensure, via bulk-boundary cor¬ 
respondence, that the boundary (namely, the edge for a 
two-dimensional system like graphene) has gapless states. 
The number of species of gapless states is given by a 
topological invariant which can be calculated from the 
bulk spectrum. While pristine graphene is gapless in the 
bulk and is therefore not topological, it can be made to 
undergo a transition to a topological phase with a non¬ 
zero Chern number by adding an appropriate conserv¬ 
ing spin-orbit coupling (SOC)^. Experimentally a SOC 


may be induced in graphene in various ways, such as by 
placing it in proximity to a three-dimensional topologi¬ 
cal insulator such as 612803^1^ or by functionalizing it 
with methyl 21 . Two models of SOC have been discussed 
in the literature: Kane-Mele^i and Rashba^^. The Kane- 
Mele type opens a gap and makes the system topological 
while the Rashba type does not open a gap and there¬ 
fore does not make it topological; consequently, in this 
work, we shall deal mostly with the former type of SOC. 
In addition, it is also interesting to consider the effects 
of an effective magnetic field with a Zeeman coupling to 
the spin of the electron. Such a coupling can arise if 
a ferromagnetic material is placed in proximity to the 
graphene^^— ; the magnetization of the ferromagnetic 
material will have only a Zeeman coupling to the electron 
spin (no orbital coupling) provided that the direction of 
the magnetization lies in the plane of the graphene. To 
the best of our knowledge, edge states, impurity effects, 
and spin transport in systems constituting spin-orbit and 
Zeeman coupled graphene have not been studied in detail 
earlier. 

In this work, we shall study the nature of edge states, 
the effects of magnetic and non-magnetic impurities, and 
spin transport in junctions involving spin-orbit coupled 
graphene both in the presence and absence of a Zeeman 
coupling term. The pristine graphene, in our work, will 
be modeled by a tight-binding lattice Hamiltonian with 
nearest neighbor hopping on a hexagonal lattice^ 

ft = - tZ E ('L(1) 

i,j OL —7,4. 

where the sum over i, j goes over the nearest neighbors, 
the hopping amplitude 7 ~ 2.8 eV, the nearest neighbor 
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spacing is d ~ 0.14 nm, and a denotes the spin compo¬ 
nent in, say, the z direction. (We will set h = \ and 7=1 
unless mentioned otherwise). The hexagonal lattice has 
unit cells which consist of two sites; we denote the upper 
and lower sites, belonging to sublattices A and B, as as 
and bfi respectively. We introduce the Pauli matrices a 
with = ±1 denoting sites on the A and B sublattices 
respectively. The midpoint of a unit cell labeled as ft is 
located at n = \/^d (ni -I- \n 2 , ^^ 2)1 where ni, 77,2 
take integer values. The spanning vectors of the lattice 
are Mi = ^d(l/2, ^3/2) and M 2 = ^d(l/2, -^3/2). 
The reciprocal lattice vectors can be chosen to be Gi = 
(47r/3d)(y3/2,l/2) and G 2 = (47r/3d)(y3/2,-1/2). As 
is well-known, such a model leads to an energy dispersion 

where 

£:- = 7|l + e*^'“i(2) 

= 7(3 -I- 2cos{V^kxd) + 4cos(^^^-^^) cos( ^^^ )]^'^^. 

The two bands touch each other at two inequivalent 
points; these are the well-known K and K' with wave 
vectors (±47r/(3-\/3d), 0). Around these points, the effec¬ 
tive low-energy continuum theory of graphene electrons 
takes the form of a (2-|-l)-dimensional Dirac Hamiltonian 
with 

Hi = - cr^fcj;)]'0g, (3) 

k 

where v = 3jd/2 is the Fermi velocity, = ±1 at K {K') 
respectively (these are called valleys), and 
denote eight-component electron annihilation operators 
with the components corresponding to sublattice (cr), val¬ 
ley (r), and spin (s) degrees of freedom. Equation ([3]) 
is the Dirac Hamiltonian and the dispersion is given by 
= ±v\k\, with a four-fold degeneracy due to the val¬ 
ley and spin degrees of freedom. 

The presence of the SOC, taken to be of the Kane-Mele 
type, and the Zeeman term arising out of proximity to a 
magnetic strip will be modeled at a lattice level by 

Hso = th ^ (c!^c.^ - cl^c.^), (4) 

Hz = 

where t 2 denotes the strength of the SOC, the sum over 
7 , j goes over next-nearest neighbors, = 1 (—1) if 
the electron makes a left (right) turn to go from site 
j to i through their common nearest neighbor, and we 
have taken the vector b = {b^, by,bz), which measures the 
strength of the effective Zeeman field, to include factors 
like the coupling to the magnetization of a proximate 
ferromagnetic material and the Bohr magneton. It is 
easy to see that Eqs. ([3]) and (O along with Eq. ([T]), lead 


to the continuum Hamiltonian near the Dirac points 

H 2 = k^-a-yky)+AsoT''a''s'"-b-^ipj:, ( 6 ) 

k 

where Ago = 3 -\/ 3 t 2 - The energy-momentum dispersion 
following from Eqs. (HD, 6 D, and ([S|) is shown in Fig. [T] 
In what follows, we shall use Eqs. (H)), (|31), and (O for 
all numerical and analytical computations done at the 
lattice level and use Eq. ([BD for analyzing the continuum 
Dirac theory for the system. 

The main results that we obtain from such an anal¬ 
ysis are the following. First, we study the edge states 
between pristine graphene (Eq. ([T|)) and graphene with 
SOC (Eq. (|3])) and demonstrate the existence of robust 
chiral edge modes provided that they are separated by 
a zigzag edge. No such modes exist for an armchair 
edge. This result is in sharp contrast to the edge modes 
between graphene with SOC and vacuum studied ear- 
lie r24,36 such modes exist both for armchair and 

zigzag edges. We also show via an exact analytical so¬ 
lution that the robustness of these edge states, in spite 
of the presence of the gapless pristine graphene, is due 
to the fact that the characteristic decay length of these 
modes vanishes in the limit O (or Ago) —>■ 0 ; this be¬ 
havior is in contrast to the usual divergence of the decay 
length edge modes with vanishing gap in the bulk. Sec¬ 
ond, we study spin-orbit coupled graphene in the pres¬ 
ence of both single and distributed impurity (impurities) 
in the weak coupling limit using a T-matrix formalism. 
We compute the energy resolved LDOS and use it to show 
that the width of the peaks in the Fourier transform of 
the LDOS provide a direct signature of the magnitude of 
the SOC. We also study a specific set of distributed im¬ 
purities and show that the corresponding LDOS reveals 
a destructive interference effect which provides a direct 
signature of the Dirac nature of graphene electrons. Fi¬ 
nally, we study the effect of magnetic impurities on the 
LDOS and show that they result in a much weaker change 
in LDOS as compared to charged impurities. Third, we 
study junctions of graphene with SOC in the form of 
both Kane-Mele (Eq. (|3])) and a Rashba term given by 
Hr = Y.ki3\[\R{T''s^(jy - sya^)]ip^, both in the pres¬ 
ence and absence oi Hz, with pristine graphene. We 
show that such junctions are necessarily spin active in 
the sense that electrons of a definite spin approaching 
a junction may reflect from it with a different direction 
of the spin. We also demonstrate that this property of 
graphene junctions may be used to generate finite, elec¬ 
trically controllable, spin currents and thus can provide 
a starting step towards applications of such junctions in 
spintronics. We note that, to the best of our knowledge, 
the presence of robust edge states, the use of LDOS in 
the presence of impurities to estimate the strength of the 
SOC, and the spin active nature of graphene junctions 
leading to finite, electrically controllable, spin currents 
in spin-orbit coupled graphene junctions have not been 
discussed in the literature. We also note that some as¬ 
pects of Kane-Mele SOC, edge states and spin transport 
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have been studied recently in buckled honeycomb sys¬ 
tems such as silicene, germanene and stanen o^^’^^ ; we 
expect our analysis demonstrating spin active junctions 
and leading to electrically controllable spin currents to 
hold for these materials as well (with minor modifica¬ 
tions to take into account the gapped Dirac spectrum of 
these materials). 

The plan of the paper is as follows. In Sec. m we dis¬ 
cuss the physics of the edge states in graphene. This is 
followed by a discussion of the LDOS due to the pres¬ 
ence of an impurity (impurities) in spin-orbit coupled 
graphene in Sec. IIIII Next, we discuss the spin active na¬ 
ture of graphene junctions in Sec. IIVI and compute the 
spin current in several possible junction geometries. Fi¬ 
nally, we discuss possible experiments, summarize our 
main results, and conclude in Sec. El 
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FIG. 1: Contour plot of energy-momentum dispersion for 
graphene with a SOC of strength t 2 = 0.05 and a Zeeman 
field bx = 0.2. E, t2 and bx are in units of 7 , while kx and ky 
are in units of 1 /d. 


II. EDGE STATES AT JUNCTION OF TWO 
REGIONS 

In this section, we study localized states at the edge 
between pristine graphene and graphene with SOC de¬ 
scribed by Eq. O and demonstrate the presence of ex¬ 
ponentially localized chiral edge states which propagate 
as plane wave along the edge. We will first consider a 
zigzag edge. 

It is known that a zigzag edge which lies between pris¬ 
tine graphene and vacuum hosts edge states for a finite 
range of the momentum along the edg o^^i^^ ; these are 


not protected by any topological symmetry since pris¬ 
tine graphene is gapless in the bulk. On the other hand, 
a zigzag edge lying between graphene with SOC and 
vacuum is known to have edge states which lie in the 
bulk gap2i; these states are topologically protected since 
graphene with SOC of the Kane-Mele type has a non-zero 
Chern number for each component of the spin. (A de¬ 
tailed discussion of edge states in graphene can be found 
in Ref. 

A system consisting of pristine graphene separated 
from graphene with SOC by a zigzag edge is gapless on 
one side. Hence if there are states on the edge, they are 
not expected to be topologically protected. However we 
will see below that for a given momentum along the edge, 
these states lie in the gap of the bulk states, which means 
they have the same momentum in both pristine graphene 
and graphene with SOC. Hence these edge states cannot 
mix with the bulk states under perturbations which con¬ 
serve the momentum. 

We are interested in studying states which are localized 
along an infinitely long zigzag edge which runs along the 
X direction. The momentum along the edge is there¬ 
fore a good quantum number. (We have absorbed the 
lattice spacing d in the definition of kx\ hence k^ is dimen¬ 
sionless). We will denote the wave functions as Am,n and 
Bm,m where the coordinate m increases vertically in the 
y direction and the coordinate n increases horizontally in 
the X direction. We assume that the wave function is a 
plane wave in the x direction so that Am,n = 
or depending on whether m is odd or 

even, and Bm,n = or depend¬ 

ing on whether m is even or odd; this is shown in Fig. [21 
We then obtain the equations (with 7=1) 


- [2 cos( )bm + bm-i] 

2^25 \s\Ti{y/^kx^o>rn sin( ^ )( — 1 

— E dm •) 

To / N 1 

[2 COS( ^ j^m E dm-\-l\ 

/ou 

+2t2S^ [sm{V3kx)bm - sm{—^){bm-i + b^+i)] 

= Ebm, (7) 


where we have taken into account the spin of the elec¬ 
tron s^. Eqs. 0 imply that we effectively have a one¬ 
dimensional system in which the site label m goes from 
—Ny/2 to Ny/2 — 1 for a finite system with 2Ny sites 
(i.e., Ny unit cells). Eqs. ([7]) will give the energy E as a 
function of the momentum kx- 

Eqs. ([7]) remain invariant under the following sets of 
transformations. 

(i) kx ^ kx + 27r/v^, a™ ->• -(-l)”^am, and bm 

{-iTb^- 

(ii) kx —>• —kx, and —>■ —s^. 

(iii) kx — >■ 27r/-\/3 — kx, E —>• —E, am —>■ (— l)”^am: and 

bm (~l)™^m- 
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Using the above transformations and combinations of 
them we can understand all the symmetries of the spec¬ 
tra shown in Figs.|3](a-d) below. [The transformations in 
(i) have a simple interpretation. The solutions of Eqs. © 
must remain invariant if the momentum is changed from 
fc to k + Gi, where Gi is one of the reciprocal lattice vec¬ 
tors given in Sec. HI Since the x component of both the 
Gi is equal to we see that Eqs. 0 must remain 

the same under ^ kx + 27r/-\/3]. 



FIG. 2: Picture of the lattice used to calculate the spectrum 
shown in Figs. [3] (a-d). The plane wave factors depend on 
the momentum k^. The label n increases horizontally in the 
positive X direction while the subscripts of am and bm increase 
vertically in the positive y direction. 


We have analytically studied these edge states lying 
between pristine graphene and graphene with SOC. We 
hnd that they take a particularly simple form if kx = 
±7r/-\/3; in this case, Eqs. o reduce to 

zb 2^25 {cLm+1 “b Clm—l') ^m —1 — ^ ^m; 
zb2t2'5 {hra-\-l “b ^m — i) = E (8) 

These equations admit a solution 

Eo = ±1/^/1+ 2\t2l 

am, bm ^ ( 1 2\t2\ ) rn>0, (9) 

where [m/2] denotes the largest integer less than or equal 
to m/2, and am = bm = 0 for m < 0. (In the limit t 2 — >■ 0, 
the wave function remains non-zero only on the four sites 
oo, bo, oi, and bi). Thus the energy lies within the bulk 
gap on the side of graphene with SOC and the wave func¬ 
tion decays exponentially on the side of graphene with 
SOC and is exactly zero on the pristine graphene side; 
this is shown in Fig. ID for t 2 = 0.1. For small t 2 , Eq. (jH]) 


We have numerically solved Eqs. 0 for a system in 
which the upper half has SOC while the lower half does 
not; more precisely, there is a SOC between sites ai and 
Qm (and between bi and bm) only if /, m are both larger 
than zero and \l — m\ < 1. The dividing line between pris¬ 
tine graphene and graphene with SOC is therefore given 
by the zigzag edge consisting of the sites Ao,„ and Bo,n 
where n goes from — oo to oo. The results are presented in 
Eigs. m (c-d) taking t2 = 0.1. For comparison, we show 
in Fig. O (a) the spectrum for pristine graphene; apart 
from the gapless bulk states which are shaded blue (they 
almost form a continuum since the bulk momentum ky 
takes a large number of almost continuous values if Ny is 
large), we see edge states which lie at exactly zero energy 
between the values of kx = —47r/(3-\/3) and —27r/(3-\/3) 
and between kx = 27r/(3-\/3) and 47r/(3-\/3). Similarly, 
Fig. [3] (b) shows the spectrum for graphene with SOC, 
taking t 2 = 0.1. The bulk spectrum is now gapped at the 
Dirac points; the gap is given by 2|Aso| = 6-\/3|t2|- We 
notice four edge states which go between the lower and 
upper bands, crossing zero energy at kx = 7r/-\/3 (see Fig. 
1 in Ref. [23) . If we now look at the spectrum shown in 
Fig. O (c) for a system with a zigzag edge lying between 
pristine graphene and graphene with SOC, we see that all 
the states present in Figs. [3] (a-b) are also present here; 
in addition, an extra set of edge states appear which lie 
very close to E = ±1 and kx = ±7r/-\/3. These are shown 
more clearly in Fig. [3] (d) which is a zoomed in view of 
the region around kx = 7r/-\/3. 


shows that the wave function decays as at a 

site which is m unit cells away from the junction inside 
the region of graphene with SOC; this implies that the 
decay length is proportional to —3(i/[ln(2|t2|/7)], where 
we have restored all the dimensionful parameters. Thus 
the decay length goes to zero as ^2 —t 0. Note that this 
behavior is in complete contrast with those of conven¬ 
tional edge states where the decay length diverges as the 
bulk energy gap vanishes. This indicates that while con¬ 
ventional edge states delocalize and merge with the con¬ 
tinuum bulk states in the limit of vanishing gap, edge 
states at the boundary of graphene with SOC and pris¬ 
tine graphene become completely localized at the zigzag 
edge separating the two regimes. Away from the special 
values of kx = ±7r/-\/3, it is difficult to obtain analytical 
solutions for the edge states. However, we find numer¬ 
ically that for small t 2 , the edge states exist only in a 
small range of values of kx close to ±7r/-\/3. The decay 
length of these states grows as kx approaches the ends 
of its allowed range beyond which the edge states merge 
with the continuum of bulk states. 

It may seem surprising that such localized edge states 
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FIG. 3: Energy-momentum dispersion for both spins (s^ = ±1) for three systems: (a) pristine graphene with a zigzag edge, 
(b) graphene with a SOC of strength t 2 = 0.1 and a zigzag edge, and (c) a system with pristine graphene and graphene with 
SOC {t 2 = 0.1) which are separated by a zigzag boundary. In all the figures, denotes the momentum in the direction along 
the edge, and the shaded regions denote the bulk states which form a continuum due to the momentum ky in the direction 
transverse to the edge. {E and k^ are shown in units of 7 and 1/d respectively). Three types of edge states are visible: (i) edge 
states of pristine graphene in Figs, (a) and (c) (these have E — 0 for a range of k^}, (ii) edge states of graphene with SOC 
in Figs, (b) and (c) (near E = 0 these have E varying linearly with \kx ± tt/v/II), and (iii) boundary states between pristine 
graphene and graphene with SOC in Fig. (c) (these lie only slightly inside the gap). Fig. (d) is a zoomed in view of the region 
around kx = tt/v^ and positive energy which shows more clearly all the edge states; the black dashed-dotted line shows the 
edge states of pristine graphene, the magenta dashed lines show the edge states of graphene with SOC, and the red dotted lines 
near the top show the boundary states between pristine graphene and graphene with SOC. 


exists even if l2 = 0 when the system has pristine 
graphene everywhere. This behavior becomes obvious 
from Eq. dH) which admits solutions with E = ±1 and 
am = ii'm-i for any value of m in this limit. Further, 
the presence of a state along a zigzag edge for any value 
of m suggests an unusually large number of states, in¬ 
creasing linearly with Ny^ at E = ±1; this is consistent 
with the Van Hove singularity in the density of states of 
pristine graphene at those two energies^. 

Finally, we have numerically studied the fate of the 


edge states for an armchair edge. We have found that 
states do not appear at an armchair edge lying between 
pristine graphene and graphene with SOC. 

We note that while graphene with SOC of the Kane- 
Mele type is a topological system with a bulk gap, pris¬ 
tine graphene is gapless in the bulk and is not a topo¬ 
logical system. Due to the gapless nature of pristine 
graphene, an edge shared between it and graphene with 
SOC is different from an edge between topologically triv¬ 
ial and non-trivial insulators. In particular, the former 
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FIG. 4: (denoting \am\'^ and \bm\'^ alternately) vs the 

coordinate m of the edge state with momentum kx = tt/V^ 
(in units of 1/d) lying at the junction between pristine 
graphene and graphene with a SOC of strength t 2 =0.1 (in 
units of 7). 


need not host localized edge states since any low-energy 
states may be delocalized since bulk pristine graphene is 
gapless. Thus one does not need to have gapless localized 
states on such a boundary; however, our work shows that 
a zigzag edge between graphene with SOC and pristine 
graphene has such states while an armchair edge does 
not. We note that these edge states do not have topolog¬ 
ical protection in contrast to their counterparts at edges 
separating topologically trivial and non-trivial insulators. 

It is interesting to compare our results on edge states 
between graphene with SOC and pristine graphene to the 
edge states which appear between (i) graphene with SOC 
and vacuum, and (ii) between pristine graphene and vac¬ 
uum. These two cases have been studied extensively in 
the literature. It is known that edge states appear at 
both zigzag and armchair edges between graphene with 
SOC and vacuum^iiSli^. This is because graphene with 
SOC is a topological system; hence states appear on any 
edge (zigzag or armchair) between this system and the 
vacuum, and all these edge states are topologically pro¬ 
tected. On the other hand, pristine graphene is a gap¬ 
less and non-topological system; hence its edges with any 
other system may or may not host any states. It turns out 
that a zigzag edge between pristine graphene and vacuum 
has edge states, but an armchair edge between pristine 
graphene and vacuum does not host any states^^— . 

We observe that Fig. 0 shows states at a zigzag edge 
for all the three cases discussed above. Figure[3](a) shows 
edge states between pristine graphene and vacuum; these 
are dispersionless and lie exactly at zero energy. Figure 


[3] (b) shows edge states between graphene with SOC and 
pristine graphene; these have a dispersion and go through 
zero energy at two particular momenta. Finally, the top 
part of Fig. [3] (d) shows edge states between graphene 
with SOC and pristine graphene; these appear only in a 
small range of momentum and have a dispersion which 
lies close to that of the bulk states. 


III. EFFECT OF LOCALIZED IMPURITY 


In this section, we will study the effect that an im¬ 
purity placed somewhere in graphene has on the LDOS 
p{f, E) as a function of the position r and energy E. By 
the LDOS we will mean the sum of the densities on the a 
and b sites at the unit cell labeled r; we will also sum over 
the electron spin. For p ristine graphene this problem was 
studied in Refs. IiTH^ our aim is to go beyond those pa¬ 
pers by studying additional characteristics in LDOS due 
to the presence of the SO term and/or Zeeman field. In 
what follows, we shall carry out an analysis of the LDOS 
in the weak impurity potential regime where perturba¬ 
tion theory holds. 

To compute the LDOS in this regime, we can use 
the standard T-matrix formalism developed for pristine 
graphene in Refs. m and HO. In the absence of any im¬ 
purities, the density of states is given by 


Po{r,E) 

Go{r,E) 


- Im {tr [Go(r,£;)]}, 

TT 


{r\ 


E - Hn 


t€ 




( 10 ) 


where the Green’s function Go(r, G) is a 4 x 4 matrix in 
sublattice and spin space, and e is an infinitesimal posi¬ 
tive number. Note that the LDOS in Eq. CHD is indepen¬ 
dent of r as a consequence of the translation symmetry of 
the system in the absence of impurities. In the presence 
of impurities, the total Hamiltonian is given by Ho + Mmp 
and the LDOS is given by 


p(r, E) 
G{f,E) 


- - Im {tr [G{r,E)]}, 

TT 

(r\ - |r). 

^ ' G - i/o - I/mp + ie ' ^ 


( 11 ) 


We now consider an impurity of strength u which is 
placed at the a site of a unit cell located at Fq with a 
potential given by 




imp 


M 4o,a afo.c- 

a=t.4 


( 12 ) 


Within the T-matrix formalism and to first order in per¬ 
turbation theory, the change in the LDOS due to the 
impurity is given by 


6p{f,E) = — — Im tr [(f^ 

TT 


1 


'^impiSo) 


E - Ho 


E — Hq -|- ie 

^ (13) 
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We define the two-point real space Green’s function for 
graphene without impurities as 


Go{ri,r2,E) 


(^il 


1 


E - Hi 
(Pk 


0 + 

t 


Ur U 

\k 


\r2) 

gifc'(ri-f2) 


( 2^)2 E - E^ 


k,cx. 


le 


(14) 


One can then write the change in the LDOS as 


Sp{r,E) = [C/o(f, tq, £1) Vimp{ro) Go{ro,r, E)]}. 

TV 

(15) 

While one can numerically compute Spif", E) using 
Eqs. (ITil) and (fTSll . the results seem to depend sensitively 
on the values of the momentum spacing l^kx, Aky and e 
that one chooses. 

In order to avoid such cutoff dependences, we have 
directly computed 6p{r, E) by numerically calculating 
p(r, E) with and without the impurity and then taking 
the difference. The calculations are carried out as follows. 
We consider a lattice in which the integers ni, n 2 go from 

I 


1 to some integer N] hence the lattice has iV^ unit cells 
and 2N'^ sites. We impose periodic boundary conditions. 
(There are two reasons for choosing such a boundary con¬ 
dition. First, it ensures that momentum is a good quan¬ 
tum number in the absence of an impurity. Second, we 
will study below the Fourier transform of the change in 
the local density of states produced by the impurity; this 
requires periodic boundary conditions in order to define 
a momentum). Corresponding to ni and n 2 , we define 
two momenta li and I 2 each of which goes from —tt to 
TT — 2'k/N in steps of 2'k/N just as we expect for orthog¬ 
onal Cartesian coordinates. In terms of the quantities rii 
and li, the phase of plane waves is given by liui + l 2 n 2 - 
We now go to the non-orthogonal position and momen¬ 
tum vectors of the hexagonal lattice by observing that 
the real space position on the lattice and the momen¬ 
tum k = (kx, ky) will satisfy hni + l 2 n 2 = kxUx + kyUy 
provided that kx = li/{'/5d) and ky = [ 2 I 2 — li)/{3d). 
Given the ranges of of li, I 2 stated above, we see that in 
the limit N ^ 00 , the Brillouin zone will be a rhombus 
with corners at (7r/d)(—1/-\/3, —1/3), {'K/d)(l/\/3,—l), 
(7r/(i)(l/-\/3,1/3), and (7r/(i)(—1/-\/3,1); the area of the 
rhombus is 87r^/(3-\/3(i^). 
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FIG. 5: Fourier transform of change in LDOS at if = 0.382 and 0.502 when an impurity is placed on the a site in one 
particular unit cell in pristine graphene (no SOC and no Zeeman field). The calculation has been done on a 30 x 30 lattice. 
(The impurity strength u = 0.1 has been divided out). The actual minimum and maximum values of the LDOS are (0,1.02) 
and (0, 0.58) respectively. The area of each picture is four times the Brillouin zone area. {E and u are in units of 7 , while kx 
and ky are in units of 1/d). 
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FIG. 6 : Fourier transform of change in LDOS at if = 0.382 and 0.517 when an impurity is placed on the a site in one particular 
unit cell, with a SOC of strength t 2 = 0.05 and a Zeeman field fea, = 0.2. The calculation has been done on a 30 x 30 lattice. 
(The impurity strength u — 0.1 has been divided out). The actual minimum and maximum values of the LDOS are (0, 1.20) 
and (0,1.58) respectively. The area of each picture is four times the Brillouin zone area. {E, t 2 , bx. and u are in units of 7 , 
while kx and ky are in units of 1 /d). 


Since we are only interested in the change in the LDOS 
to first order in the impurity strength u, we will take 
M to be a small number and calculate Sp{r,E)/u. We 
first consider pristine graphene (i.e., in the absence of 
SO coupling and Zeeman field) when an impurity of 
strength u = 0.1 is placed at the a site of the unit 
cell centered at fp = (0,0). The absolute value of 
the Fourier transform of Sp(f’,E)/u is shown in Fig. [S] 
for different values of E, with a 30 x 30 lattice (i.e., 
N = 30). Since the energy spectrum is found to have 
an exact or nearly exact six-fold degeneracy at most 
energies, we will calculate the LDOS by summing over 
the contributions from the six states with energy clos¬ 
est to the desired value of E. Figure [S] shows that we 
get large and sharp peaks at the Dirac points, namely, 
the six points forming a hexagon around the center. 


[The six points are given by (47r/3-\/3(i) times (1,0), 
(—1/2, •\/3/2) and (—1/2, —•\/§/2), and (47r/3-\/3d) times 
(—1,0), (1/2, \/3/2) and (1/2, —-^3/2) which are respec¬ 
tively equal to K and K' up to reciprocal lattice vectors]. 
The peaks broaden as we move away from zero energy. 
(We note that the normalization of the LDOS calculated 
in this way is arbitrary to the extent that we have done 
the calculations for a particular system size and have not 
normalized the results to take that into account). We 
also note that the Fourier transform of Sp{f, E) is always 
zero at fc = (0, 0) since that is just the difference in the 
number of states at that energy with and without the 
impurity, and we have chosen E in such a way that the 
impurity does not change that number. 

In Fig. ini we show the absolute value of the Fourier 
transform of 6p{r,E)/u when an impurity of strength 
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(a) (b) (c) 


FIG. 7: Fourier transform of change in LDOS at E = 0.382 and 0.517 when an impurity of strength 0.1/6 is placed at each of 
the six sites around one particular hexagon (hence the integrated strength is equal to 0.1). There is a SOC of strength t 2 = 0.05 
and a Zeeman field = 0.2. The calculation has been done on a 30 x 30 lattice. (The impurity strength u = 0.1 has been 
divided out). The actual minimum and maximum values of the LDOS are (0, 0.03) and (0,0.05) respectively. The area of each 
picture is four times the Brillouin zone area. [E, t 2 , bx and u are in units of 7 , while kx and ky are in units of 1/d). 



FIG. 8 : Fourier transform of change in LDOS at FI = 0.331 and 0.483 when a magnetic impurity (coupling to = ±1 with 
strengths ±0.1 respectively) is placed on the a site in one particular unit cell, with a SOC of strength t 2 = 0.05 and no Zeeman 
field bx- The calculation has been done on a 30 x 30 lattice. (The impurity strength u = 0.1 has been divided out). The actual 
minimum and maximum values of the LDOS are (0, 0.004) and (0, 0.007). The area of each picture is four times the Brillouin 
zone area. {E, t 2 , bx and u are in units of 7 , while kx and ky are in units of 1/d). 


M = 0.1 is placed at the a site of the unit cell centered at 
r = (0, 0), when there is a SOC of strength t 2 = 0.05 and 
a Zeeman field = 0.2. Comparing Figs. [5] and [SI we 
find that the SOC and Zeeman field broaden the peaks 
at the Dirac points. This is expected since the SOC and 


Zeeman field open a gap and broaden the Dirac points 
in the dispersion shown in Fig. [T] Thus we find that 
the breadth of the LDOS peaks at the Dirac point is a 
measure of the strength of the SOC and/or Zeeman field 
in graphene. 
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Next, we study a distribution of non-interacting im¬ 
purities instead of a single fixed impurity studied ear¬ 
lier. More specifically, we place impurities of strength 
Ui = 0 . 1/6 on each of the six sites around a graphene 
hexagon, so that the integrated impurity strength is 
u = X]i=i = 0.1 as before. We compute the absolute 
value of the Fourier transform of 6p{r,E)/u, as shown 
in Fig. [71 for t 2 = 0.05 and a Zeeman field = 0.2. 
Comparing Figs. [5] and [71 we see that there no peaks at 
the Dirac points when the SOC is present. Further, the 
maximum value of the Fourier transform is now much 
smaller than in the case when the impurity is present 
only on a single site. These observations can be un¬ 
derstood as follows. The Fourier transform of the two- 
point Green’s function in Eq. o is particularly large 
at the Dirac momenta K and K' . Equation (fT51) shows 
that Sp{f,E) is composed of two such Green’s functions. 
Hence the Fourier transform of Sp{f, E) will be peaked at 
the difference of K and K' , namely, at K and K' (since 
K — K' = K'), provided that the Fourier transform of 
Vimp does not vanish at K and K' . This is true if there is 
an impurity at a single site. However, if there are impuri¬ 
ties of equal strengths at the six sites around a hexagon, 
the Fourier transform of this, given by ’'S van¬ 

ishes at k = K and K' due to destructive interference 
between contribution from each point. Hence the Fourier 
transform of 5p{r, E) is negligible at K and K—-^. Such 
a cancellation is unique for Dirac electrons in graphene 
and has been pointed out in the context of LDOSi^ and 
STM spectral^ of a single impurity placed at the hexagon 
center in graphene; our work here points out that such 
a cancellation is qualitatively important for understand¬ 
ing the structure of LDOS for distributed impurities in 
graphene. 

Finally, we study the LDOS in the presence of a mag¬ 
netic impurity at a single site of graphene with SOC. 
Such an impurity provides a direct coupling to electron 
spin at that site. More specifically, we assume that the 
impurity is on the a site of a unit cell located at tq and 
couples with strength u to the x-component of the spin: 

Vimp = U (at^ -|- Q'f'pqflro.t)- (16) 

For u = 0.1 and a SOC of strength t 2 = 0.05, the Fourier 
transform of the change in the LDOS is shown in Fig. |5| 
for two values of the energy E. Comparing Figs. [5] and 
|8l we Hnd that the scale of the change in the LDOS 
is much smaller for a magnetic impurity compared to 
a non-magnetic impurity of the same strength, namely, 
0.004—0.007 versus 1.20—1.58. This can be partly under¬ 
stood as follows. Since and s“ anticommute, a unitary 
transformation of the Hamiltonian by leaves the SOC 
parameter t 2 unchanged but flips the impurity parameter 
u ^ —u. Since the LDOS must be invariant under this 
unitary transformation, it must be an even function of 
u. To lowest order, therefore, the change in the LDOS 
must be of order for a magnetic impurity, while it is 
of order u for a non-magnetic impurity. For u = 0.1, we 


therefore expect the change in the LDOS to be about 10 
times smaller for a magnetic impurity. Thus we find that 
a magnetic impurity will have a smaller effect on LDOS 
compared to its non-magnetic counterpart. 


IV. SPIN ACTIVE GRAPHENE JUNCTIONS 

In this section we will study the differential conduc¬ 
tance G for either a junction of graphene with SOC 
and pristine graphene or two regions of pristine graphene 
which are separated from each other by finite width strips 
of various kinds, such as graphene with SOC or in an 
external Zeeman field. In Sec. IIV A1 we carry out an an¬ 
alytical calculation for the differential conductance from 
a continuum theory. This will be followed by Sec. IIVBI 
where we will numerically calculate G for finite-sized sys¬ 
tems using a lattice model. A comparison between the 
results obtained by these two approaches is given in Sec. 

wu 


A. Analytical calculation using continuum models 



FIG. 9: Schematic representation of the junction geometry 
studied in Sec. IIV Al The shaded region represents graphene 
with SOC while the white region denotes pristine graphene. 

In this section, we analyze transport in graphene junc¬ 
tions with SOC. The geometry for such junctions which 
will be studied in this section is shown in Fig. |9l We be¬ 
gin with an analysis of the geometry in Fig. |9| (a) which 
represents a junction of pristine graphene and a part of 
graphene which has both Kane-Mele and Rashba SOC 
terms. 

To analyze transport across such a junction, we first 
consider the system shown in Fig. |5| (a); the junction lies 
at y = 0. In region H where y > 0, the Hamiltonian is 
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given by 

k 

= V {T^a^k^ + la^dy) + Aso 

+ Xr {T^a^sy - (jys^), (17) 

where the momentum kx is a conserved quantity hav¬ 
ing the same value everywhere. In the presence of the 
Rashba term, the energy-momentum dispersion is given 
by a quartic equation for E, 

[E^ - v\kl + kl) - /Xlf = {E-A,,f. (18) 

The solution of Eq. m leads to a gapped energy spec¬ 
trum with four energy bands as shown in Fig. [TUlfor rep¬ 
resentative values V = 3/2, t 2 = 0.05, Xr = Aso/10, and 
ky = 0. We observe that the spectrum is not symmetric 
about E = 0. 



FIG. 10: Energy vs kx in a region of graphene with Kane- 
Mele SOC of strength t 2 = 0.05, Rashba SOC of strength 
Xr = Aso/10, and ky = 0. {E, t 2 and Aso are in units of 7 , 
while kx is in units of 7 and 1 /d). 


Given some values of E and kx (which remain the same 
in the regions of pristine graphene and the strip region 
with SOC), the momentum ky in the strip can take four 
values given by 

K = - Al, ± 2Xr\E - Asol, (19) 

where the ± sign outside the square root is independent 
of the ± sign inside. We thus have four possible values of 
the momentum ky. Depending on the different parame¬ 
ters some of these values may be imaginary. If they are 
imaginary we will consider only the exponentially decay¬ 
ing solutions, while if they are real, we will choose the 
signs so that the group velocity dE/dky is positive so 
that the electrons are moving right, i.e., towards y = 00 . 


In any case, only two out of the four possible values of 
ky are physically allowed in region II; let us denote these 
two values by ky’^. 

In what follows, we further use the fact that is a 
good quantum number. We will therefore only study 
the case = 1. The case = — 1 gives similar re¬ 
sults since it is related to = I by the unitary trans¬ 
formation hj: —>■ T^ayhp-^(jy. The operator is not 
a good quantum number. However, we observe that 
h{—kx) = h{kx)<jy. Since this transformation flips 
both and kx , it is enough to study the case of an inci¬ 
dent electron with = I and all values of kx- 

In region I of Fig. IH] (a), the Hamiltonian is given by 
Eq. (ITTI) with Ago = A_r = 0. In this region, is also a 
good quantum number. The wave functions for right and 
left moving spin-up (s^ = 1) £md spin-down (s^ = —I) 
electrons with momentum {kx,iky) (where ky > 0) and 

energy E = vJkx + ky are given by 


^i{±kyy+kxX — Et) 


i{±kyy+kxX — Et) 
^ 1 


where e®“ = {kx — iky)/^kx + ky. We now consider 
a spin-up electron which is incident on the junction 
with momentum {kx,ky) and energy E = fi, where 
/i = Er + eV is the chemical potential or voltage applied 
in region I measured with respect to the Dirac point and 
Er is the Fermi energy. The reflected wave function can 
then be written as Tpr = where r^-t' 

rj,-]- are functions of E and kx- Note that represents 
the amplitude for an incident spin-up electron to be re¬ 
flected from the junction as a spin-down electron. It is 
therefore a direct measure of the spin active nature of 
the junction. Such a reflection process which converts a 
spin-up electron to a spin-down electron constitutes an 
analog in spin space of Andreev reflection from a super¬ 
conductor in which an incident electron is converted to 
a reflected hole. The total wave function in region I can 
thus be written as 

■0/ = V’-i-r + (21) 

In region H, the presence of the Rashba term implies 
that is not a good quantum number. Consequently, 
the transmitted wave function will have amplitudes in 
both = 1 and = — 1 sectors. For transmitted elec¬ 
trons with energy E and momentum {kx, ky) or {kx, ky), 
the wave functions can be found by solving the equation 
h^|;a = E^a, where a = 1,2. Note that ky'"^ can be real 


/ 1 


■^±1 = - 7 = 




V2 0 

V 0 


V'it = 


1 


/ 0 

0 
1 
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or imaginary. A straightforward calculation yields 




UAt 

Mat 

Uai 


m 


(\ 

,.a 

V J 

ia°‘ 


pi{Ky+kxX-Et) 


2Xr{E-Aso) el’ 
ia°' 


2XRel’ 
E - A,,„ 


‘‘Bi 


= 1 , 


- v\kl + kl) - AL, 

v{ka;-ik^), el = v{k^ + iky), 


( 22 ) 

(23) 


current Xf and the reflected spin-down current are 
given by 



We show plots of TZ\,^ as a function of the applied voltage 
^ for a fixed Xr = A^,o/10 in Fig.[TT](a), and as a function 
of Xr for a fixed /i = 1 in Fig. [TT] (b). These clearly 
demonstrate the spin active nature of the junction. We 
see that for a fixed Xr, TZfi indeed vanishes at /r = A^o, 
but it eventually increases with /i. This demonstrates 
that the spin current can be electrically controlled. 


where is a normalization constant which ensures that 
tpl^tljl = 1. (The value of is not required in the ex¬ 
pressions presented below). The transmitted wave func¬ 
tion in region II is thus given by tpu = X)a=i 2 ^aV'+- 
To find r-i-i-, and we impose continuity of the 
wave function at the junction: 'tpi{y = 0 ) = ipiiiy = 
0). This leads to the following conditions on the various 
amplitudes: 


1 -f r^'l' 

71 

72 

72 


— 'y ^ taU^A, 
0 = 1,2 

= y ^ 

0 = 1,2 

“ y ^aUAli 
0 = 1,2 

= ^ ^aUBl- 

0 = 1,2 


(24) 


The solution to these equations yields 


V 


(l-e- 2 .) (j^_A,,) (e^-e 7 

e]_el V ’ ^ ^ 

+ {u\^u\^ + “Ai^Bt ~ ~ “Ai^Bt)- (2®) 


Note that a non-zero value of is a consequence of the 
presence of two solutions, with ky = ky’^, for a fixed en¬ 
ergy E and transverse momentum kx', these two solutions 
merge when A/j = 0 and the junction ceases to be spin 
active in this limit. We observe that r 4 ^'t' = 0 for all val¬ 
ues of kx if E = Ago; this shows that the specific voltage 
at which the spin-flip transport takes place can be con¬ 
trolled by the gap originating from the Kane-Mele term. 
(rj^-|- = 0 also vanishes if the incident electron comes in 
at a glancing angle, namely, if = 0 so that e*“ = ± 1 ). 
For the incident electron, the range of values of kx goes 
from — fco to ko, where fco = since we want Ek,,,Ay 

to be equal to with real values of ky. Integrating over 
this range of kx, we find that the total incoming spin-up 
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FIG. 11: Reflection probability TZi^ from a junction between 
pristine graphene and graphene with a Kane-Mele SOC of 
strength t 2 = 0.05 and a Rashba SOC of strength Xr. Figure 

(a) shows TZi'f as a function of /r for Xr — Aso/10, while figure 

(b) shows 7?.^^ as a function of Ah for r = 1. {r and Xr are 
shown in units of 7 ). 
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Next, we analyze the geometry shown in Fig. |9] (b). 
In this geometry, pristine graphene resides in regions I 
and III, and graphene with Kane-Mele and Rashba SOC 
forms an interface region II of width d which lies between 
those two regions. To analyze the transport in this sys¬ 
tem, we note that in regions I and III, is a good quan¬ 
tum number. In region I, the wave function is given by 
Eq. (I2ID . Similarly, in region III, the transmitted wave 
function is given by 


i’ll! = 


(28) 


where denotes the probability for a spin-a electron to 
be transmitted when the incident electron has = 1. 

In region II, the electron wave function is a linear su¬ 
perposition of electrons with momenta The wave 

function of an electron with momentum ky’’^ is given 
by Eq. while that for an electron with momentum 
— ky''^ is given by 


V'- = 

= 

v%^ = 

= 


1 

w 


( \ 


Vb^ 


\^bJ 


^i{-k^y-\-kxX-Et) 


la 


2Xr{E-Aso) 


2A«e“ 
E - 


Vbi 


= 1 , 


(29) 


where is a normalization constant which ensures 
= 1- Using Eqs. (0^ and (1^^ . the wave func¬ 
tion in region II can be written as 


i’ll = X] + dai’-)- (30) 

a^l,2 


yields 

a^l,2 

iPaV-B^ + QaVB-x)^ 

a^l,2 

a^l,2 

X + gqUgj^), (31) 

a=l,2 

a^l,2 

a^l,2 

a^l,2 

a^l,2 

(32) 

Using these we can compute the reflection and transmis¬ 
sion probabilities Ra^ = \ra^\'^ and Ta-f = |tatP respec¬ 
tively, where a can be f or J,. Similarly, if the electron 
incident from region I had = —1, we would have the 
reflection and transmission probabilities Rai and T^j,. 
These must satisfy the unitarity relations 

T.|-|- -I- -I- -I- = 1, 

U ~ 1- (33) 

In what follows, we shall compute the tunneling con¬ 
ductance by solving Eqs. (13111321) . and compute the trans¬ 
mission probabilities between two pristine graphene re¬ 
gions (region I and III) across a strip of graphene (re¬ 
gion II) with a width do = 15d, for t 2 = 0.05 and 
Xr = Ago/lO. In Figs. [T^] (a) and (b), we show the 
transmission probabilities and 

versus kx for two values of the energy E. Our plots 
clearly demonstrates a finite spin conversion as indicated 
by the dashed blue lines in Figs. m a) and (b). Given 
the transmission probabilities, the differential conduc¬ 
tances can be calculated as follows. For a momentum 
k = {kx,ky), the current in the y direction is given by 
Jy ^ = \dEi^/dky\. Let the chemical potentials in regions 
I and III in Fig. [H](b) be and /i 2 , so that the voltage 
bias between the two regions is given by eV = {y. 2 —pi)/ e. 
In the zero bias limit in which F —> 0 and /ii, ^2 P, 
the differential conductance Ga^ = dl/dV for an inci¬ 
dent electron with spin a being transmitted with spin /3 
(a, (3 can be f or 1) is given by 


1 -I- 

~vr 

e“ -f 

72 

72 

72 

~7^ 

72 

~7^ 

^^^gi(fc„d+a) 

71 


The amplitudes Pa, qa, io-t: and ra-^ can be found by 
matching the wave functions at a; = 0 and x = d. This 


G^p{p) = e^IU 


dkxdky 


d{p-Ej:)T^0j^j:, (34) 
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where W is the width of the system in the x direction (we 
assume that W ^ d). Integrating the (5-function over ky 
in Eq. (IMl) gives a denominator equal to \[dE^/dky) 
which precisely cancels the ^ appearing in the numer¬ 
ator of that equation. We thus obtain 

g2w rko 

GaflifJ-) = J7rS2 / ( 3 ^) 

J — ko 

where ko = \n\/v. 

Instead of plotting Gaisifk) versus /i, it is convenient to 
plot the ratio Gajs (fk) /Go (m) ; where Go (/i) is the conduc¬ 


tance when there is perfect transmission, i.e., Tap = Sap- 
From Eq. (1351) we find that Go(/i) = [e^VE/(27r)^][2|/r|/?;] 
in one particular valley. We then have the expressions 


G-|''|' 

Go 


1 

2ko 

1 

2fco 



dkx T-f-f, 
dkx 


(36) 


Fig. [12] (c) shows plots of G-|-^/Go and G^^/Go versus fi. 



(a) 



(b) 



FIG. 12: Transmission probabilities and conductance across a strip of graphene with a width of 15 (in units of the lattice 
spacing d), a Kane-Mele SOC of strength t 2 = 0.05, and a Rashba SOC of strength Xr = Aso/10. (a) and (b) show the 
transmission probabilities T-f^ (solid red) and (dashed blue) vs for E — 0.4 and 0.5 respectively, while (c) shows G^^/Go 
(solid red) and (dashed blue) vs p. {t 2 and Aso are in units of 7 , while kx is in units of 1/d). 


Another interesting quantity to consider is the rota¬ 
tion of the electron spin produced by region II. For each 
value of E and kx, we know that a spin-up electron in¬ 
cident from region I converts to a linear superposition of 
spin-up and spin-down on being transmitted to region III, 
with amplitudes and respectively. In spin space, 
the linear superposition (tf-i-, t 4 ,t)^ describes an electron 
whose spin polarization points at an angle 6 with respect 
to the z axis, where tan(0/2) = We can there¬ 

fore define an average rotation angle produced by region 
II as 

(0) = / dkx 0 (% + T^t): (37) 

^^0 J-ko 

where we have weighted the angle of rotation by the 
transmission probability -b T|-|-. In Fig. |T3|we show 
the average rotation angle as a function of p for transmis¬ 
sion across a strip of graphene with the same parameters 


as in Fig. [T^l Fig. |l3| clearly shows a finite spin rota¬ 
tion which increases as a function of p in the zero-bias 
limit; this demonstrates the potential of these junctions 
as generators of electrically controllable spin current. 

Before ending this section, we note that the calcula¬ 
tions of this section serve as a proof of principle that a 
finite electrically controllable spin current may be gen¬ 
erated in graphene junctions with finite SO coupling. 
There are many ways to enhance the magnitude of this 
current, which in our chosen parameter regime, appears 
to be rather small. For example, one can increase the bias 
voltage V and/or the thickness d of region II which will 
increase both and 9. Further, there may be other, 
more suitable geometries for larger spin current gener¬ 
ation. In the next sub-section, we shall carry out nu¬ 
merical calculation from a lattice Hamiltonian which will 
address some of these issues and verify the approximate 
continuum calculation of the present section. 
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FIG. 13: Average rotation angle vs /r for transmission across 
a strip of graphene with a width of 15 (in units of the lattice 
spacing d), a Kane-Mele SOC of strength t 2 = 0.05, and a 
Rashba SOC of strength Xr — Aao/10. {fj., t 2 and Aso are in 
units of 7 ). 


B. Numerical calculation using lattice models 

In this section, we provide a numerical analysis of the 
transport across the various junctions starting from a lat¬ 
tice model. The advantages of a lattice calculation over 
a continuum calculation are that a lattice calculation is 
straightforward to implement numerically (for instance, 
one does not have to impose any matching conditions on 
the wave functions), and one can easily study the effects 
of potentials or spin-orbit couplings which have arbitrary 
spatial profiles. 

We will follow a procedure similar to the study of con¬ 
ductance across junctions of graphene and other materi- 
al o^'^^ . We will assume that the strip is infinitely long in 
the X direction, so that the momentum of an electron 
incident from one of the regions of pristine graphene is 
a good quantum number everywhere in the system. The 
incident energy E is also a good quantum number. How¬ 
ever, ky will vary from one region to another depending 
on the presence of SOC and a Zeeman field. 

We will calculate the conductance numerically using a 
lattice model similar to the one shown in Fig. [51 We con¬ 
sider an electron incident from the pristine graphene at 
the bottom of that figure and we calculate the probabil¬ 
ities of reflection (back to the bottom) and transmission 
(to the pristine graphene at the top). An incident spin-up 
electron can get either transmitted or reflected as spin-up 
or spin-down; we will denote the corresponding probabili¬ 
ties by T^-f, and respectively. Similarly an 
incident spin-down electron will have transmission and 
reflection probabilities given by T 74 ,, R^i, and R^. 

The calculation is done as follows. Given the values of 
the momentum k^ and energy E (which we will hence¬ 
forth assume to be positive), the dispersion for pristine 

I 


graphene given in Eq. m uniquely fixes a momentum ky 
lying in the range [0, 27r/3d]. (It may happen that there 
is no real solution for ky] this would imply that such a 
value of E is not allowed for the given momentum k^. In 
that case we will set the transmission probabilities equal 
to zero). Then the incident and transmitted waves will 
have momentum ky while the reflected wave will have 
momentum —ky. We now consider a single transmitted 
wave, with unit amplitude and equal to either 1 or — 1 , 
which is located at the top of Fig. [51 and we find which 
superposition of the four possible incident and reflected 
waves at the bottom would give rise to such a transmit¬ 
ted wave (we have to allow for four possible waves in 
general since they could be either incident or reflected 
and they could have = ±1). This superposition can 
be found by using Eqs. © to set up a matrix problem 
where the four reflection and incident amplitudes as well 
as the values of am and bm inside the region with SOC 
or Zeeman field appear on the left side of an equation 
and the single transmitted wave at the top (with unit 
amplitude) appears as a source term on the right of the 
equation; the reflection and incident amplitudes are then 
found by doing a matrix inversion. Having found these 
amplitudes for the two cases where the transmitted wave 
has equal to 1 and — 1 , we then invert these relations 
and find the reflection and transmission amplitudes when 
a wave is incident with unit amplitude = ±1. The 
modulus squared of the amplitudes give the reflection 
and transmitted probabilities as usual. Einally we check 
if the unitarity relations in Eq. ( 1551 ) are satisfied. 

Given the transmission probabilities, the differential 
conductances can be calculated as described in Sec. lIV Al 
We again arrive at Eqs. ((34l) and (l35l) . except that the 
range of integration of kx in the lattice model is given by 
[—27r/-\/3, 27r/-\/3]. However only those values of kx will 
contribute for which E^ can be equal to /r with real values 
of ky. Once again, we will plot the ratio Gai 3 {fJ-)/Go{^)^ 
where Go(/i) is the conductance when = Sa/ 3 . Given 
a chemical potential ^ lying between 0 and 1 , we can 
show using Eq. that 


Go(m) 


e^IF 


J dkx 

8 /1-Mn 


acos( 



(38) 


If /i is small, Go{fj.) varies linearly with /i, namely, 
Go(/r) = [e^kF/(27r)^][8/x/u] where v = 3/2 is the Fermi 
velocity. This expression is exactly twice of what we ex¬ 
pect for two species (due to the valleys) of massless Dirac 
electrons in the continuum. The additional factor of two 
is because we have considered the full range of kx from 
[—27r/-\/3 to 27r/\/3]; this double counts the contribution 
from each of the two valleys since the transmission is in¬ 
variant under kx ^ kx + 211 j^/i. The double counting is 
not present in the ratio Gap{^)/G q{^). 
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FIG. 14: Transmission probability T-t-t- vs kx across a strip of graphene with 15 unit cells and a SOC of strength t 2 — 0.05. The 
values of energy are (a) E = 0.3, (b) E = 0.4, and (c) E = 0.5. {E and t 2 are in units of 7 , while kx is in units of 1/d). 


We now present our numerical results for a number 
of different cases. We first consider a strip of graphene 
with a width of 15 unit cells and a SOC of strength t 2 = 
0.05; there is pristine graphene on both sides of the strip. 
(We take the Zeeman held to be zero). We will study 
the transmission probabilities Tap as a function of the 
incident energy E and the momentum kx- Since the SOC 
does not couple spin-up and spin-down electrons, we will 
have = T-fp = 0. Further, using the symmetries 
discussed after Eqs. o, we can show that 

Tnikx) = T^-fikx+ 2Tr/y^), 

Ttt(fc,) = Tui-kx). (39) 

The second equation in (15^ implies that it is sufficient 
to study In Fig. 1141 we show the transmission prob¬ 
ability as a function of kx for three values of the 
energy, E = 0.3, 0.4 and 0.5. In each of the hgures we 
see that there are regions of kx where Tfl exactly zero 
or close to zero. These regions occur for two reasons. 
First, we have already seen that in pristine graphene, 
for a given value of kx, all possible values of E are not 
allowed; for a disallowed value of E, we set r|-|- = 0. 
Second, in graphene with a SOC of strength t 2 = 0.05, 
the minimum value of energy occurs at the four values 
kx = ±27r/(3-\/3) and ±47r/(3-\/3) and that minimum en- 
ergy is given by Ago = 3v^t2 — 0.26. In Fig. [14] (a), 
the energy E = 0.3 is only a little bit more than Ago. 
Hence for all values of kx except the regions around the 
four special momenta, the energy of the incident elec¬ 
tron lies inside the gap of graphene with SOC, and the 
wave function will decay exponentially inside that part 
of graphene. is therefore very small for all values of 
kx except near those four momenta. 

In Fig. m we see some transmission resonances, 
namely, for certain values of E and kx, we find that 



FIG. 15: G^t/Go vs p for transmission across a strip of 
graphene with 15 unit cells and a SOC of strength t 2 — 0.05. 
(/r and t 2 are in units of 7 ). 


is close to 1. We can understand this as fol¬ 
lows. If the energy does not lie in the gap of region 
of graphene with SOC, i.e., if the momentum ky = 

±{l/vWE^-v^kl-Al in that region is real, then we 
expect a transmission resonance if (3/2)kyNy is an inte¬ 
ger multiple of tt. This is because such a condition implies 
that the wave function in the region with SOC will sat¬ 
isfy tp{y = Ny) = ±tp{y = 0), where the ± sign depends 
on whether {3/2)kyNy is an even or odd multiple of tt. 
Hence the wave function will match at j/ = 0 and Ny 
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between pristine graphene and graphene with SOC, with 
the reflection amplitude being equal to zero at ?/ = 0 and 
the transmission amplitude being equal to ±1 at y = Ny] 
we will therefore get = 1. 

In Fig. [151 we show G^^/Gq as a function of /r. We see 
that G^t/Go is very small for ji < Ago — 0.26. As ^ is in¬ 
creased to 1, G^-t-ZGo also approaches 1 although some os¬ 
cillations are visible. The locations of the maxima can be 
qualitatively understood as follows. We saw in the pre¬ 
vious paragraph that there are transmission resonances 
if {2>/2)k'yNy = UTT, where n = 1,2, 3, • • •. Since Fig. [T3] 
shows that the resonances are most prominent close to 
the Dirac points kx = ±27r/(3-\/3) and ±47r/(3-\/3), let 
us ignore the contributions from values of kx away from 
the Dirac points and approximate the dispersion inside 

graphene with SOC hy E = ^J{vky)^ + this holds if 

E is not too large. We therefore expect Gf-|-/Go to show 
maxima when 

The smallest values of n = 1,2,3 give ^ = 

0.334,0.493,0.680 which are approximately the locations 
of the first three maxima in Fig. [121 

The system discussed above, with only a SOC present, 
enjoys an additional symmetry, namely, 

Tnik.) = Tuikx). (41) 

In fact, even the transmission amplitudes are equal, 
= tii{kx)- This can be shown as follows. We first 
note that for a particular value of equal to either 1 or 
— I, Eqs. (O have a symmetry resembling time reversal 
in which all numbers are complex conjugated. (This does 
not change the value of kx which simply appears as a pa¬ 
rameter in those equations. This symmetry is therefore 
a bit different from the usual time reversal symmetry in 
which both kx and ky change sign). For a particular value 
of equal to 1 or —1, this implies that the scattering 
matrix S which relates the incoming waves at the top and 
bottom of the system to the outgoing waves must be sym¬ 
metric, in addition to being unitary. (This can be proved 
as follows. If ii and i 2 denote the incoming amplitudes at 
the top and bottom, with plane wave factors and 

gi 3 fc„/ 2 ^ and oi and 02 denote the outgoing amplitudes at 
the top and bottom, with plane wave factors and 

e“*3fcy/2^ they must be related as { 01 , 02 )'^ = ■ 

Complex conjugating this relation transforms 1*^2 to O1/2 
and vice versa. Time reversal symmetry then implies that 
we must have = 'S'(oi,02)^- This implies that 

S't = S'*, namely, S is symmetric. Hence the transmis¬ 
sion amplitude S'21 from the top to the bottom must be 
equal to the transmission amplitude S12 from the bottom 
to the top. Next, we use the fact that the Hamiltonian 
in Eq. ([HI) (but without a magnetic field b) is symmetric 
under the parity transformation y ^ —y, x ^ x, namely, 

^ a" s". (42) 


The transformation in Eq. (I42p interchanges the a and b 
sublattices and also flips the component of the spin. 
This symmetry implies that the transmission amplitude 
from the top to the bottom for an electron with spin 
must be equal to the transmission amplitude from the 
bottom to the top for an electron with spin —s^. Com¬ 
bining these two symmetries, we see that the transmis¬ 
sion from the bottom to the top must be the same for 
= ± 1 . 



FIG. 16: r_>_» (solid red) and T<_>/Go (dashed blue) for 

transmission of an electron with energy E = 0.5 and = 1 
across a strip of graphene with 15 unit cells and a SOC of 
strength t2 — 0.05; the bottom seven unit cells have V = —0.2 
and the top eight unit cells have V = 0.2. {E, t 2 and V are 
in units of 7 , while kx is in units of 1/d). 

If the parity symmetry is broken, by applying a poten¬ 
tial Vm which depends on the ?/—coordinate m in a way 
which is not invariant under reflection about the center of 
the region with SOC {Vm can be independent of both the 
spin and the sublattice index), we expect that the trans¬ 
mission amplitudes and will no longer be equal. 
Taking linear combinations of the incident electron so as 
to be quantized along, say, the x direction, we And that 
the transmission amplitudes for a = I electron to be 
transmitted as a = I and —1 electron are given by 
{t^^ + ti^)/2 and (t-|-|- — tii)/2 respectively. (In fact these 
expressions hold for any component of the spin which is 
perpendicular to the z axis, not just s^). The latter will 
not be zero in general which implies that the s“ com¬ 
ponent of the electron can flip when it transmits across 
a region with SOC and a parity-breaking potential. We 
demonstrate this effect in Eigs.fTOlandlTTlfor transmission 
across a strip of graphene with a width of 15 unit cells 
with a SOC of strength t 2 = 0.05; in addition, the bot¬ 
tom seven unit cells have a potential V = —0.2 and the 
top eight unit cells have V = 0.2. Eig.lTSl shows the prob¬ 
abilities for an incident electron with energy E = 0.5 and 


hkx, — k.. 
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FIG. 17: G->->/Go (solid red) and G<—>/Go (dashed blue) 
vs n for transmission across a strip of graphene with 15 unit 
cells and a SOC of strength t 2 = 0.05; the bottom seven unit 
cells have V = —0.2 and the top eight unit cells have V = 0.2. 
(/i, t2 and V are in units of 7, while kx is in units of 1 /d). 

= 1 to be transmitted into an electron with = ±1; 

the transmission probabilities and _are shown 

as functions of kx- (Here —denote = ±1). We 
see that there is a non-zero (though small) probability of 
conversion from s“ = 1 to —1. Figure ITTI shows plots of 

and _^ versus /i for the same system. (The spin 

conversion effect discussed here is related to spin filter 
and spin valve effects which have been discussed in other 
papers, for instance, Refs. [29I - I3TI) . 

Next we consider a strip of graphene with a width of 
15 unit cells and a Zeeman field in the x direction with 
bx = 0.1. (We take the SOC strength to be zero). To 
study this problem, we generalize Eqs. © to include a 
Zeeman field bx which couples spins f £md 4-5 now 

Til, ^nd Til ^.ll be non-zero in general. We can 
then derive some symmetries similar to the ones discussed 
after Eqs. ©; using these we find that 

Tafiikx) = Ta0{kx +2Tr/\/3), 

Tafjikx) = T^a-pi-kx), (43) 

where we define —a =4, (f) if a =t (4-) and similarly 
for —/3 in the second equation. It is therefore enough to 
study r.|-|- and Tn. (If t 2 = 0, we also have the symmetry 
Tapikx) = Tap{-kx)). 

In Eig. [181 we show Gh/Gq and Gh/Gq as functions 
of fj,. We see that Gh/Gq is much larger than G-t-|-/Go 
in the entire range of fi. This happens for this partic¬ 
ular value of Ny = 15 and can be qualitatively under¬ 
stood as follows. If electrons with = ±1 (rather than 

= ±1) were incident, they would be transmitted with 
unit magnitude but with a phase difference. If the elec¬ 
trons have energy E, the dispersion inside graphene with 



FIG. 18: Gii/Go (solid red) and Gj_-|-/Go (dashed blue) vs n 
for transmission across a strip of graphene with 15 unit cells 
and a Zeeman held = 0.1. (p and are in units of 7). 


a Zeeman field would be given by if = (vky±) ± bx, if 
E is not too large. (For simplicity, we are assuming 
that kx is equal to one of the Dirac points ±27r/(3-\/3) 
or ±47r/(3-\/3)). Here ky± denote the values of the y 
component of the momentum for = ±1. The phase 
difference between the electrons with = ±1 is given 
by (3/2)(fcj^_ — ky+)Ny = 3bxNy/v. For bx = 0.1 and 
Ny = 15, the phase difference is 3 which is close to tt. 
Hence electrons with = ±1 are perfectly transmitted 
but with almost opposite signs. Hence incident electrons 
with = 1, which is given by the linear combination 
(|s“ = 1) -|- |s“ = —1))/-\/2 will be transmitted almost as 
the linear combination (js^^ = 1) — |s^ = —1))/-\/2 which 
is the same as = —1. We thus see an almost perfect 
conversion of spin from = 1 to —1. Note that this 
approximate argument is independent of the energy E 
which explains why Gh/Gq is much larger than G^^/Gq 
for all y in Fig. [181 

Finally we consider a strip of graphene with a width 
of 15 unit cells and a SOC of strength = 0.05, fol¬ 
lowed immediately by another strip with the same width 
of 15 unit cells where there is a Zeeman field bx = 0.1 
but no SOC; there is pristine graphene on both sides of 
the two strips. We consider an incident electron which 
first strikes the region with SOC and we study the trans¬ 
mission after it leaves the region with a Zeeman field. 
Once again T77, Tn, Tn, and Tn will all be non-zero in 
general, and the symmetries in Eqs. (H51) will hold. In 
Fig. [m we show Gh/Gq and Gh/Gq as functions of y.. 
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FIG. 19: G^^/Go (solid red) and Grt/Go (dashed blue) vs n 
for transmission across two successive strips of graphene, the 
first one with 15 unit cells and a SOC of strength t 2 = 0.05, 
and the second one with 15 unit cells and a Zeeman field 
63 , = 0 . 1 . (/r, t 2 and are all in units of 7 ). 


C. Discussion of spin active jnnctions 

To summarize and compare the results presented in 
Secs. II V Al and lTV^ we have discussed two kinds of junc¬ 
tions which are spin active, i.e., they can rotate the di¬ 
rection of spin of an electron which is incident on the 
junction. The first example, discussed in Sec. IIV Al is 
a region of graphene which has a combination of Kane- 
Mele and Rashba SOC. The Rashba SOC does not con¬ 
serve the spin; hence it is not unexpected that it can 
give rise to a spin active junction of spin-orbit coupled 
and pristine graphene. We have used a continuum the¬ 
ory (valid near the K and K' points) to analytically cal¬ 
culate the reflection probability from a junction of spin- 
orbit coupled and pristine graphene and the transmission 
probability and differential conductance (obtained by in¬ 
tegrating the transmission over all incident momenta) 
through a strip of spin-orbit coupled graphene. To quan¬ 
tify the spin active nature, we have studied the amount 
of spin rotation as a function of the applied voltage and 
the strength of the Rashba SOC. In the second example, 
discussed in Sec. IIVBI we studied the effect of a strip of 
graphene with Kane-Mele SOC, a Zeeman field (in a di¬ 
rection perpendicular to the SOC so that the two terms 
do not commute), and a potential which is not parity 
symmetric. Since the calculation cannot be analytically 
done for a general non-parity symmetric potential, we 
have used the tight-binding model to numerically calcu¬ 
late the transmission probability and difference conduc¬ 
tance across such junctions. Once again we find that the 
junction is generally spin active. 

Comparing the results for the different junctions, we 


see that a Zeeman field which is perpendicular to the 
Kane-Mele SOC and a Rashba SOC are most effective in 
producing spin active junctions. A non-parity symmetric 
potential along with a Kane-Mele SOC is relatively less 
effective. 


V. DISCUSSION 

In this work, we have studied edge states, effects of im¬ 
purities, and spin active junctions in spin-orbit coupled 
graphene along with the presence/absence of a Zeeman 
term which originates from proximity of the graphene 
sheet to a suitably chosen ferromagnetic film. The SOC, 
which may arise due to proximity of graphene to topolog¬ 
ical insulator films, has been taken to be of either Kane- 
Mele or Rashba form. 

Our study concentrates on three properties of such 
graphene systems. First, we have shown that a junction 
between graphene with the Kane-Mele type of SOC and 
pristine gapless graphene (with no SOC) supports robust 
chiral edge states provided that the edge separating the 
two regions is of the zigzag type; no such states exist 
for the armchair edges. We have also shown that these 
edge states are robust in spite of the presence of gapless 
pristine graphene on one side of the junction; this robust¬ 
ness arises due to the fact that the decay length of these 
states vanishes for t2 —>■ 0. We have pointed out that 
such behavior is in complete contrast to the behavior of 
conventional edge modes where the decay length diverges 
in the limit of vanishing gap. 

Second, we have studied the change in the LDOS orig¬ 
inating from either a single or a specific distribution of 
impurity atoms in spin-orbit coupled graphene. We have 
shown that for a single impurity the Fourier transform 
of the LDOS displays peaks near the Dirac points with 
a finite width; the width of these peaks is a direct mea¬ 
sure of the strength of the induced spin-orbit interaction. 
We have also shown that for a specific distribution of im¬ 
purity atoms (distributed at the corners of a graphene 
hexagon), the Fourier transform of LDOS exhibits an ab¬ 
sence of peaks near the Dirac points. Such an absence 
can be traced back to the destructive interference of the 
contribution to the LDOS from each of the impurity sites 
and is a direct signature of the Dirac nature of graphene 
electrons. Such an effect has been discussed earlier in 
the context of LDOS^^ and STM spectra of single im¬ 
purity placed at the center of a hexagon in graphene^S; 
however, its manifestation has not been pointed out for 
a distribution of impurities to the best of our knowledge. 

Although we have only discussed the effects of a single 
impurity or a small number of impurities in this paper, 
our results can also be used to understand what would 
happen if there was a finite density of impurities which 
are far from each other, so that the scattering from the 
different impurities is incoherent. The Fourier transform 
of the change in the LDOS would then be given by the 
Fourier transform of the LDOS for a single impurity mul- 
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tiplied by the density of impurities. Hence the Fourier 
transform of the LDOS of a finite density of impurities 
will share the features of the Fourier transform of the 
LDOS of a single impurity such as the peaks at the Dirac 
points. The Fourier transform of the LDOS for a finite 
density of impurities can be measured by a light scatter¬ 
ing experiment. 

Third, we have studied junctions of spin-orbit cou¬ 
pled graphene (with both Kane-Mele and Rashba terms) 
and pristine graphene. We have shown that such junc¬ 
tions are generally spin active and that they may be 
used to generate electrically controllable spin currents 
in graphene. We have demonstrated this in a variety of 
junctions with analytic computations using low-energy 
effective Dirac-like Hamiltonians and with numerical cal¬ 
culations based on microscopic lattice models. We have 
also discussed several ways of enhancing the spin current 
and pointed out the role of Zeeman coupling terms and 
parity-symmetry breaking potential terms in this con¬ 
text. 

The experimental verification of our work would in¬ 
volve preparation of graphene samples with strong SOC. 
Since the intrinsic SOC in graphene is extremely weak, 
this needs to be done using a proximate material with 
strong SOC; the hybrid samples of topological insulators 
atop a graphene sheet which have already been experi¬ 
mentally studied are ideal for this purpose. The LDOS 
for impurities in such samples can be measured using 
an STM; the Fourier transform of the LDOS can then 
be computed^. The prediction of our present work is 
that the width of the peak of the Fourier transform of 
this LDOS would be a direct measure of the strength of 
the induced SOC. To form junctions of spin-orbit coupled 
graphene with its pristine counterpart, we need to deposit 


the topological insulator over a part of the graphene sam¬ 
ple leaving the rest of the sample in its pristine form. We 
predict that if such a junction has a zigzag edge separat¬ 
ing the spin-orbit coupled and pristine graphene, there 
would be additional chiral edge states whose density of 
states could be measured by STM; no such edge states 
would exist if an armchair edge separates the two regions. 
Finally, for spin active junctions we suggest measurement 
of the spin current via standard tunneling conductance 
measurements where the injection and detection of cur¬ 
rent is done with spin-polarized leads of opposite polari¬ 
ties. This will serve as a direct measure of G-j-j,. 

To conclude, we have studied edge states, effects of im¬ 
purities, and spin active junctions in graphene which has 
spin-orbit coupling. Our results points out the presence 
of robust chiral edge states in graphene junctions sepa¬ 
rating spin-orbit coupled graphene from its pristine coun¬ 
terpart with novel properties of their decay length, shows 
that the local density of states originating from impuri¬ 
ties in spin-orbit coupled graphene near the Dirac points 
can serve as a measure of the strength of the induced 
spin-orbit coupling for graphene, and demonstrates that 
junctions of spin-orbit coupled and pristine graphene are 
spin active and may be used to generate electrically con¬ 
trollable neutral spin currents. We have proposed realis¬ 
tic experiments which may test our theory. 
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